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Section  1 

INTRODUCTION 


Side-chain  thermotropic  liquid  crystalline  (LC)  materials  have  been  investigated  for 
possible  applications  in  electro-optical  devices1,  nonlinear  optics2-3,  selective  optical  filters  and 
reflectors4-5,  and  optical  storage  media5-7.  In  particular,  linear  polysiloxane  polymers  have 
received  considerable  attention,  as  a  promising  backbone  moiety  due  to  their  low  glass  transition 
temperatures.  The  flexibility  of  this  backbone  minimizes  the  dynamic  coupling  of  the  main-  and 
side-chain  motions.  Recently,  several  groups  have  also  reported  ov.  possible  microphase 
separation  of  the  hydrocarbon  side-chain  mesogens  from  the  linear  siloxane  backbone  based  on 
experimental  X-ray  diffraction  results8-10.  Electron  density  calculations8  indicate  dial  the  smectic 
layer  can  be  divided  into  sublayers  which  consist  of  a  mesogenic  core,  spacer  groups,  and  aliphatic 
tails  with  the  siloxane  backbone  forming  a  thin  (5  A-7  A)  interface  between  layers. 

Liquid  crystalline  systems  of  increasing  interest  consist  of  small  siloxane  rings  with 
different  pendant  mesogenic  groups  attached,  for  example,  a  mixed  cyanoester  and  cyanobiphenyl 
system11.  Another  such  LC  is  based  on  cyclic  penta(methylsiloxane)  with  combinations  of 
cholesteryl-4'-allyloxybenzoate  {C}  (40%  mole  fraction)  and  biphenyl-4'allyloxybenzoate  (B) 
(60%  mole  fraction)  mesogens  (cf.  Figure  1).  The  packing  and  structure  exhibited  by  these 
commercially  available  cholesteric  LC's  (Wacker  LC-Silicones)12  have  been  recently  studied, 
showing  unusually  well-defined  interlayer  order  giving  rise  to  high  order  X-ray  reflections13.  This 
LC  is  of  technological  importance  because  of  attributes  related  to  light  processing  media5-6. 

Molecular  mechanics  and  dynamics  are  well  established  techniques  for  raacromolecular 
structural  definition14  at  the  atomic  level,  applied,  for  example,  for  polymeric  chains153’6.  A 
molecular  simulation  of  a  cyclic  siloxane-based  liquid  crystal  model  by  employing  semiempirical 
and  molecular  mechanics  calculations  was  recently  performed16.  This  computational  study  offered 
an  insight  into  the  effects  of  molecular  structure  on  the  intramolecular  interactions  and  also  the 
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Figure  1:  "Hie  chemical  structure  of  the  studied  cyclic  siloxine  -  based  liquid  crystal. 
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Section  2 

THEORY 


2.1.  MOLECULAR  MODELING 

The  molecular  mechanics  and  molecular  dynamics  modeling  methods  are  employed  here  to 
study  the  packing  geometries  of  a  cyclic  siloxane-based  LC.  In  order  to  perform  a  molecular 
dynamics  simulation,  one  must  first  describe  tbe  potential  energy  surface  of  the  molecular  system 
in  question.  The  potential  energy  surface  is  obtained  in  this  case  by  the  application  of  the  semi- 
classical  molecular  mechanics  (MM)  approach.  An  MM  calculation  of  a  molecule  evaluates  its 
conformational  energy,  and  an  energy  minimization  represents  it  as  it  would  exist  near  0°  K.  By 
simulating  a  molecule  at  elevated  temperatures  for  the  calculation  of  structural  fluctuations,  phase 
transitions  and  time  averaged  properties,  the  molecular  dynamics  application  is  used  to  partition  the 
total  energy  between  kinetic  and  potential  contributions,  and  to  distribute  the  latter  between  the 
vibrational  states  of  tbe  molecule.  Molecular  mechanics  and  molecular  dynamics  calculations  in 
this  study  were  performed  by  using  the  Quanta/CHARMm17  molecular  modeling  computer 
software  on  a  Silicon  Graphics  workstation. 


2.1.1.  Molecular  Mechanics 

2.1.1.1  Potential  Energy 

In  the  molecular  mechanics  methodology  the  potential  energy  of  a  molecule  or  group  of 
molecules  is  expressed  in  a  classical  approach  as  the  sum  of  terms  describing  the  state  of  the 
system.  These  include  bond  lengths,  bond  angles,  torsion  angles  and  nonbonded  pairwise 
interactions.  The  potential  energy  of  a  system  of  configuration  r  may  then  be  written  as  follows: 
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V(r)=  l/2£Kb(b-b0)2+  1/2ZK0  (8-e0)2+l/2£K<p(l4cos  (n<p-«))+£(A/r12-C/i6+qlq2/Dr)  (1 ) 

bonds  bond  angles  torsion  angles  non-bonded  pairs 

where  b,  0,  <p,  and  r  denote  bond  lengths,  bond  angles,  torsion  angles  and  nonbonded 
intramolecular  distances,  respectively.  The  force  field  parameters  Kb,  Kfi,  K(p,  A,  C,  and  D,  the 
geometrical  values  bo,  ©o,  n,  8  and  the  partial  atomic  charges  q  are  described  for  the  CHARMm15 
force-field  in  Reference  18.  The  parametrization  for  the  organosiloxane  derivatives  are  taken  from 
Reference  19,  and  the  partial  atomic  charges  for  the  siloxane  ring  from  Reference  16. 

In  order  to  calculate  the  potential  energy  of  a  molecule,  or  group  of  molecules,  by  the  use 
of  equation  (1),  a  set  of  (x,  y,  z)  cartesian  coordinates  is  required  to  describe  the  position  r  of  each 
atom.  For  example,  the  contribution  to  the  potential  energy  expression  from  the  bond  lengths 
would  be  calculated  as  described  in  the  following.  Fust,  the  bond  length  bij  between  two  bonded 
atoms  is  calculated: 

bij  =  {  (xi  -  xj)2  +  (y,  -  yj)2  +  (z;  -  zj)2  }  M*  (2) 

where  {xi,  yj,  z;},  and  {xj,  yj,  zj),  are  the  cartesian  coordinates  for  atoms  i  and  j,  respectively. 
The  standard  (optimum)  value  for  this  particular  bond  length  (bo)  is  then  subtracted  from  the  value 
determined  in  equation  (2).  In  order  to  obtain  the  potential  energy  of  this  bond  length,  the  (b  -  bo)2 
term  is  multiplied  by  an  appropriate  force  constant,  Kb-  The  procedure  is  repeated  for  all  of  the 
bonds  present  within  the  molecule17.  Such  calculations  are  performed  for  all  energy  terms  in 
equation  (1).  A  decrease  in  computation  time  of  the  nonbonded  atom  pair  interactions  term  is 
obtained  by  introducing  a  cut  off  distance  (8  A  - 15  A  depending  upon  the  molecular  size). 

2.1. 1.2  Energy  Minimization 

Once  the  potential  energy  surface  of  the  molecule  is  described,  its  conformation  can  be 
varied  in  an  effort  to  reduce  the  potential  energy  of  the  system  and  better  describe  the  lowest  energy 
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geometry.  Several  energy  minimization  techniques  can  be  invoked  ior  energy  minimization.  The 
simplest  procedure  (and  often  most  useful)  is  the  Steepest  Descent  method  which  involves  the 
calculation  of  the  gradient  of  the  potential  energy  of  the  system  with  respect  to  the  coordinates. 
Specifically,  an  iterative  process  is  performed  in  which  the  gradient  of  the  potential  is  calculated 
with  respect  to  a  particular  set  of  coordinates.  A  new  value  of  the  coordinates  is  then  found  which 
lowers  the  potential  energy  of  the  system.  A  constant  (typically  1.2)  is  then  multiplied  by  the 
negative  gradient  of  the  potential  and  a  new  energy  calculated.  If  the  energy  decreases,  the  process 
continues.  If  the  energy  increases  after  a  step  has  been  taken,  the  step  size  is  reduced  by  a  factor 
(typically  0.5)  for  the  next  iteration20.  Equation  (3)  describes  one  iteration  (i)  of  this  procedure  for 
the  minimization  of  a  torsion  angle  q>: 

(pi  =  (pi-f-ki dV/ap  (3) 

Although  the  Steepest  Descent  method  is  suitable  for  removing  major  distortions  in 
geometry  or  improper  steric  interactions,  it  may  locate  local  potential  energy  minima  (rather  than  a 
global  minimum)  for  the  molecular  geometry21.  Thus,  more  powerful  energy  minimization 
techniques  have  been  developed  such  as  the  Conjugate  Gradient  method  and  the  Newton-Raphson 
method,  discussed  extensively  elsewhere21. 

2.1.2.  Molecular  Dynamics 

2.1. 2.1  General 

Molecular  dynamics  calculations  simulate  the  natural  motion  of  a  molecular  system  and 
produce  a  set  of  trajectories  which  describe  the  atomic  motions  of  the  system  over  time21. 
Calculations  can  be  performed  at  constant  temperature  or  constant  pressure.  They  yield  various 
time  averaged  properties  of  the  system,  such  as  potential,  kinetic  or  total  energies  of  the  system, 
torsion  angles,  interatomic  distances  and  bond  angles.  In  addition,  the  trajectories  show  the 
“instantaneous”  values  of  these  quantities  as  a  function  of  time.  In  general,  molecular  dynamics 


calculations  have  proved  to  be  useful  in  examining  the  behavior  of  molecular  systems  over  time  at 
specified  temperatures  or  pressures. 


2. 1.2.2  Methodology 

Consider  a  system  of  N  atoms,  where  atom  i  has  mass  m;  and  position  rj.  The  relationship 
between  the  velocity  (dri  /  dt)  and  the  momentum  (pj)  of  the  atom  is  given  by: 

drj/dt  =  pi/mi  (4) 

Substituting  this  into  Newton's  equation  of  motion,  F  =  ma,  yields: 

Fi  =  dpi  /  dt  (5) 

where  F;  is  the  net  force  exerted  on  atom  i  by  all  of  the  other  atoms  in  the  system.  Since  the  force 
is  the  negative  of  the  gradient  of  the  potential  energy,  equation  (5)  can  be  rewritten  as: 

Fi  =  -3V/9rj  =  dpi  /  dt  (6) 

The  force,  F;  is  obtained  from  the  gradient  of  die  potential  energy  (equation  (1)).  B>  integrating 
equation  (6)  over  small  enough  time  steps  (At  -  1/2  fs),  the  velocity  and  position  of  each  atom  can 
be  updated  to  a  new  velocity  and  position  by  using  the  Veriet  numerical  integration  method22.  The 
procedure  yields  a  trajectory  for  the  molecular  system  over  some  specified  time  period. 

The  molecular  dynamics  calculation  may  be  performed  at  constant  temperature  by  making 
use  of  the  following  relationship: 

3  kb  T  =  X  mi  vj.v;  /  N  (7) 
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where  kb  is  the  Boltzmann  constant  In  order  to  keep  the  temperature  of  the  system  constant,  the 
velocities  of  the  atoms  are  adjusted  appropriately  (scaled)  using  a  method  developed  by 
Anderson23 . 

It  is  appropriate  to  note  that  in  most  cases,  molecular  dynamics  calculations  require  a  large 
amount  of  computer  time.  In  order  to  reduce  computation  time  the  SHAKE  algorithm25  may  be 
employed,  which  allows  the  user  to  increase  the  time  step  used  in  the  calculation  ai  j  thereby 
reduce  the  total  number  of  iterations  to  complete  one  simulation.  The  time  step  used  must  be  small 
compared  to  the  period  of  the  highest  frequency  motion  of  interest  in  the  system.  Because  C-H 
stretching  is  typically  the  highest  frequency  of  vibration  in  a  molecule,  SHAKE  is  used  to  fix  these 
bonds  at  their  original  equilibrium  values.  So  doing,  the  time  step  used  in  a  molecular  dynamics 
simulation  can  be  increased  by  a  factor  of  2,  i.e.,to  a  value  of  1  femtosecond  (fs). 

2.1.2.3  Simulation  Procedure 

A  molecular  dynamics  simulation  consists  of  heating,  equilibration  and  simulation  steps, 
briefly  described  below. 

2.I.2.3.I.  Heating 

The  first  step  in  a  molecular  dynamics  simulation  is  obtaining  a  set  of  cartesian  coordinates 
for  the  atoms  constituting  the  molecule  or  group  of  molecules.  From  these  the  potential  energy 
surface  is  specified  by  using  the  MM  potential  function  described  in  equation  (1).  In  addition, 
prior  to  the  initiation  of  a  molecular  dynamics  calculation,  the  potential  energy  of  the  molecular 
system  should  be  partially  minimized  to  remove  major  distortions  in  geometry  or  improper  steric 
interactions. 


This  minimized  starting  structure,  with  all  atoms  having  zero  velocities,  represents  the 
molecular  system  at  a  temperature  near  absolute  zero.  In  order  to  perfoim  a  molecular  dynamics 
simulation  at  a  higher  temperature,  it  is  necessary  to  heat  the  molecule.  Heating  is  generally 
accomplished  by  slowly  and  uniformly  adding  kinetic  energy  to  the  molecular  system  at  periodic 
time  intervals  by  assigning  velocities  to  each  atom,  either  randomly  or  in  proportion  to  the  force 
acting  on  the  atom.  The  length  of  the  heating  period  is  extremely  important  and  is  dependent  upon 
the  size  of  the  molecular  system.  If  the  heating  period  is  too  short,  achieving  a  constant 
temperature  may  prove  to  be  impossible  in  the  equilibration  and  simulation  stages.  For  a  system 
containing  approximately  700  atoms  (such  as  some  of  those  studied  here),  a  heating  time  of 
approximately  2  pico  seconds  is  sufficient.  Smaller  systems  can  be  heated  faster  and  larger 
systems  should  be  heated  longer. 


2.1.2. 3.2.  Equilibration 

After  the  system  is  heated  to  the  desired  temperature,  it  is  necessary  to  equilibrate  it  over  a 
period  of  time  sufficient  to  make  the  total  energy  of  the  system  (potential  +  kinetic)  remain 
constant.  This  procedure  uses  as  input  the  dynamics  files  created  from  the  heating  phase  of  the 
simulation.  Equilibration  is  continued  until  the  total  energy  becomes  constant  over  time. 


2.I.2.3.3.  Simulation 

After  the  system  has  been  equilibrated  at  the  desired  temperature,  the  simulation  phase  is 
carried  out  for  10  ps  or  more  in  order  to  get  reasonable  time  averaged  results  and  other  statistical  or 
thermodynamic  information.  The  simulation  step  uses  as  input  the  dynamics  files  created  during 
the  equilibration  phase  of  the  calculation. 
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2.2.  X-RAY  DIFFRACTION  CALCULATION 

Cylindrically  averaged  scattering  is  the  most  common  technique  for  modeling  scattering 
from  an  anisotropic  material,  namely  oriented  with  (cylindrical)  symmetry,  without  assuming  any 
crystalline  order2®.  Averaging  is  performed  by  rotation  around  the  axis  of  the  direction  of  the 
fibrous  microscopic  model  (by  convention  the  z  axis),  so  that  the  intensity  transform  I(ft,  Z)  is 
cylindrically  symmetrical  about  it  This  simulation  is  most  applicable  for  highly  oriented  materials, 
such  as  in  the  LC  state2®-27.  The  Fourier-integral  of  the  scattering  amplitude  for  n  atoms  is : 

n 

F(S)  =  Xfj(s)e2ltl(S-rj)  (8) 

j 

where  fj  indicates  the  scattering  power,  S  the  reciprocal  lattice  vector,  and  rj  the  interatomic  vector 
in  real  space.  It  can  be  shown  that  this  can  be  expressed  in  terms  of  the  n-th  order  Bessel  function 
Jn  for  a  solution  in  cylindrical  coordinates^ : 

l(RZ)  =  Z  (£  fj  J„  (ftrj)  cos(Zzj  -  n<bj))2  +  £  (1  fj  Jn  (ftrj)  sin(Zzj  -  nfy))2  (9) 

n  j  n  j 

where  (ft  ,  Z,  y;  X=ftcosy,  Y=ftsiny,  Z=Z  )  and  (r,  z,  x=rcos$,  y=rsin(>,  z=z),  are  the 
cylindrical  coordinates  in  reciprocal  and  direct  space,  respectively,  and  I  =  IF12.  Note  that  I  is  only 
a  function  of  ft  and  Z  for  the  cylindrical  symmetrical  component,  being  averaged  with  respect  to  y 
for  the  oriented  sample.  The  summation  in  equation  (9)  is  done  for  all  atoms  j,  and  -5<n<+5, 
since  only  these  terms  of  the  Bessel  function  are  of  importance2®-27. 

Spurious  features  which  may  occur  in  this  general  diffraction  simulation  of  molecular 
models  are  eliminated  by  using  a  model-size  correction;  in  particular,  the  cylindrically  averaged 
scattering  is  subtracted  from  that  derived  for  a  prism,  or  cylinder,  with  the  same  radial  structure  as 
the  model  but  structureless  in  the  z  -  direction.  Such  a  calculation  describes  a  special  case  of  the 
above  general  formulation  in  equation  (9)2®.  This  total  scattering,  including  the  model-size 
correction  as  well  as  the  meridional  section  (namely  the  Fourier  transform  of  the  projection  of  the 
structure  onto  the  z  axis)  were  calculated  by  the  CERIUS  simulation  program27. 
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Section  3 

RESULTS  AND  DISCUSSION 


3.1.  Dynamics  on  Isolated  Molecules 

An  understanding  of  the  conformational  flexibility  and  degree  of  interdigitation  in  the  cyclic 
siloxane-based  LC  shown  in  Figure  1,  may  be  initially  obtained  by  molecular  dynamics 
calculations  on  the  isolated  molecular  models.  The  three  models  are  those  in  which  all  mesogens 
are  either  axial  or  equatorial  (shown  in  Figures  2  and  3,  respectively),  or  at  an  intermediate 
conformation  (shown  in  Figure  4)  in  which  two  {C}  sidechains  are  in  the  axial  position,  and  three 
{ B )  sidechains  are  in  the  equatorial  position.  In  these  three  figures,  the  sizes  of  the  silicon  atoms 
are  not  proportional  to  those  of  the  other  atoms.  These  have  been  termed  the  cone,  disk  and 
cylinder  models13,  respectively,  because  of  their  characteristic  geometries,  and  have  been 
previously  studied16.  Details  of  the  molecular  dynamics  simulations  are  presented  in  Table  la. 

Structural  differences  between  the  starting  models  (Figures  2-4)  and  those  from  the 
molecular  dynamics  calculation  (Figures  5-7)  were  examined.  The  structures  shown  in  Figures  5- 
7  are  the  lowest  energy  and  averaged  structures  of  the  simulation  after  a  10  pico  second  ran  at  300° 
K.  Clearly,  a  comparison  of  Figures  2-4  and  5-7  shows  a  great  deal  of  conformational  variations. 
Further  analysis  was  carried  out  by  comparing  various  torsional  angles  during  the  simulation  to 
those  in  the  initial  structure. 

The  mean  variation  of  a  torsion  angle  during  the  molecular  dynamics  run  is  derived  from  a 
histogram  plot  of  the  dihedral  angle  distribution  during  the  simulation,  from  which  the  mode  can  be 
determined.  This,  in  turn,  approximates  the  mean  since  a  normal  distribution  is  generally 
observed.  As  an  example.  Figure  8  shows  a  histogram  for  one  of  the  five  dihedral  angles  within  a 
mesogen.  Specifically,  Tori ,  the  first  rotatable  bond  of  the  mesogen  {O-Si-Cfint  of  spacer-mesogen- 
Csecond  of  spacer-mesogen)  is  plotted.  Tables  II  and  III  summarize  the  mean  magnitudes  and  range  of 
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TABLE  I  :  Details  of  Molecular  Dynamics  Simulation3 


Cal  Isolated  Cone.  Disk  and  Cylinder  Models 


Coneb 

Diskb 

Cylinder0 

Heating  time  (ps) 

0.5 

0.5 

2.5 

Equilibration  time  (ps) 

64.0 

5.5 

41.5 

Simulation  time(ps) 

10.0 

10.0 

10.0 

Time  step  (ps) 

0.0005 

0.0005 

0.001 

Equilibration  time  (ps)  1 30.0 

Simulation  time(ps)  10.0 

Time  step  (ps)  .  0.0005 


Heating  time  (ps) 

- 515 - 

Equilibration  time  (ps) 

40.0 

Simulation  time(ps) 

10.0 

Time  step  (ps) 

0.001 

Heating  time  (ps) 

5.0 

Equilibration  time  Ops) 

10.0 

Simulation  time(ps) 

10.0 

Time  step  (ps) 

0.001 

3  Non  bonded  cut  off  distance :  8  A;  non  bonded  pair  list  update  frequency :  every  50  iterations. 
b  Hie  SHAKE  procedure  was  not  applied. 
c  The  SHAKE  procedure  was  applied 
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figure  5:  Model  structures  of  the  cows  isomer  after  a  10  picosecond  molecular  dynamics 
simulation.  The  upper  structure  (a)  is  the  lowest  energy  conformation  adopted  during  the 
simulation.  The  lower  structure  (b)  is  the  statistically  averaged  structure  of  the  molecule  during  the 

10  pico  second  simulation. 
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(*) 


(b) 

Figure  6:  Model  structures  of  the  disk  isomer  after  a  10  pico  second  molecular  dynamics 
simulation.  The  upper  structure  (a)  is  the  lowest  energy  conformation  adopted  during  the 
simulation.  The  lower  structure  (b)  is  the  statistically  averaged  structure  of  the  molecule  during  the 

10  pico  second  simulation. 
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Figure  7:  Model  structures  of  the  cylinder  isomer  after  a  10  pico  second  molecular  dynamics 
simulation.  The  upper  structure  (a)  is  die  lowest  energy  conformation  adopted  during  the 
simulation.  Tbe  lower  structure  (b)  is  the  statistically  averaged  structure  of  the  molecule  during  the 

10  pico  second  simulation. 
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the  Tori  angles  in  the  three  models  from  the  10  pico  second  dynamics  simulation.  Table  II  mostly 
relates  the  initial  model  to  the  “equilibrated”  structure,  while  Table  HI  relates  to  deviations  from  this 
equilibrated  structure  during  the  MD  run. 


TABLE  II  :  <Torl> *  and  <ATorl>b  Mean  Dihedral  Angles  (°)  from  the  Molecular 


Mesogen 

Type 

Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

<Torl> 

<ATorl> 

<Torl> 

<  ATorl  > 

<Torl> 

<ATorl> 

(C> 

Initial 

60 

63 

185 

MD 

55 

5 

74 

11 

176 

9 

{C} 

Initial 

45 

59 

300 

MD 

68 

23 

54 

5 

283 

17 

{  B  } 

Initial 

180 

54 

60 

MD 

153 

27 

65 

11 

66 

6 

(B) 

Initial 

180 

175 

61 

MD 

284 

104 

195 

20 

64 

3 

{B  } 

Initial 

40 

55 

69 

MD 

187 

147 

71 

16 

67 

2 

aTorl  -  {O-Si-Cfustofspacer-mesogen-Csecondofspacer-mesogen)- 
b  The  mean  variation  of  the  torsion  as  compared  to  the  initial  value. 


TABLE  III  :  Range  of  the  Tori  Dihedral  Angles  (°)  during  the  Molecular 
Dynamics  (MD)  Simulation 


Mesogen  Type 

Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

{C> 

49 

88 

56 

(C) 

50 

62 

38 

(B> 

50 

170 

84 

(B) 

56 

76 

58 

(B) 

42 

74 

41 
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No  of  Structures 


Notably,  the  largest  conformational  flexibility  around  this  bond  is  observed  in  the  cone 
model,  with  two  of  the  dihedral  angles  (on  {B}  mesogens)  varying  by  104°,  and  147°  on  average 
from  the  initial  magnitudes,  while  all  other  <Torl>  differences  are  smaller  than  30‘.  The  large 
variations  of  this  dihedral  angle  may  reflect,  in  part,  the  stability  of  a  particular  model  with  respect 
to  the  proximity  of  the  mesogens.  Thus,  the  result  is  anticipated  since  the  cone  model  with  all  of 
the  mesogens  at  an  axial  position  is  the  least  stable  energetically  (the  minimum  energies  of  the 
initial  cone,  disk,  and  cylinder  models  are  309,  298,  and  260  kcal/mol,  respectively)16.  The 
average  <ATorl>  values  are  61°,  13°,  and  7°  for  the  cone,  disk,  and  cylinder  models,  respectively. 
This  relative  flexibility  of  the  cone  model  is  also  reflected  in  the  range  values  shown  in  Table  Ill. 
The  values  for  the  cone  model  are  on  average  (=  49°)  much  smaller  than  those  for  the  disk  model 
(=  94°),  indicating  a  lower  conformational  flexibility  of  the  mesogens.  This  may  be  due  to  an 
increase  in  the  steric  interactions  among  the  mesogens  for  the  cone  relative  to  the  disk  model. 
Interestingly,  although  the  cylinder  model  has  the  lowest  minimum  energy  of  the  three16,  its 
mesogens  have  on  average  a  lower  conformational  flexibility  than  those  in  the  disk,  as  also 
evidenced  by  a  lower  average  range  of  torsion  angles  (=  55°). 

A  different  pattern  of  variations  in  these  models  is  observed  for  the  Tort  torsions  (defined 

by  the  position  of  the  mesogen  relative  to  the  ring  :  |  Si-O-Si-Cfirsi  of  spaccr-mesogcn))  which 

describe  possible  flexibility  in  the  orientation  of  the  mesogens  relative  to  the  siloxane  ring 

Statistical  data  of  the  results,  tabulated  for  the  cone,  disk  and  cylinder  models  in  Tables  IV  and  V, 

indicate  the  largest  conformational  flexibility  around  this  bond  is  observed  fot  the  cone  model  with 

all  of  the  five  dihedral  angle  differences  in  the  range  of  ±120*.  The  average  \<ATor2>\  values  arc 

121°,  75°,  and  71°  for  the  cone,  disk,  and  cylinder  models,  respectively.  Specifically,  two  of  the 

mesogens  which  had  an  initial  Tort  torsion  of  -60*  with  respect  to  the  siloxane  ring  assume  mean 

values  of  ca.  180*,  thus  resembling  the  cylinder  model.  Once  again,  this  may  indicate  a  lower 

stability  for  the  cone  model.  In  general,  however,  all  of  the  three  models  exhibit  large  changes  in 

the  mesogen's  orientation  relative  to  the  siloxane  ring.  The  range  of  torsion  angles  shown  in 

Figures  9(a)-(c)  indicate  that  both  the  cone  and  cylinder  model  ranges  are  approximately  60°  for  all 

21 


angles,  while  a  larger  range  is  observed  for  the  disk  model.  These  changes  can  be  clearly  seen  in 
Figure  6.  The  greatest  possibility  for  conformational  flexibility  in  the  position  of  the  mesogens 
relative  to  the  siloxane  ring  is  exhibited  by  the  disk  model.  This  can  be  explained,  in  pan,  by 
lower  steric  interactions  between  mesogens  in  this  case,  enabling  easy  reorientation. 

TABLE  IV  :  <Tor2> 4  and  <ATor2>b  Mean  Dihedral  Angles  (°)  from  the  Molecular 
_ Dynamics  (MD)  Simulation  in  Comparison  to  the  Initial  Values 

Mesogen  Isolated  Cone  Isolated  Disk  Isolated  Cylinder 

Type 


<Tor2> 

<£Tor2> 

<Tor2> 

1 

V 

<Tor2> 

<ATor2> 

{C} 

Initial 

304 

176 

184 

MD 

184 

120 

295 

119 

189 

5 

{C} 

Initial 

306 

176 

309 

MD 

172 

134 

49 

127 

59 

1 10° 

{B} 

Initial 

303 

194 

298 

MD 

63 

120° 

183 

11 

169 

129 

{B> 

Initial 

305 

176 

293 

MD 

59 

114c 

176 

0 

282 

11 

{B} 

Initial 

299 

54 

175 

MD 

181 

118 

293 

121° 

77 

98 

aTor2-  (Si-O-Si-Cfirst of spacer-mcsogcn ) ■ 

b  The  mean  variation  of  the  torsion  as  compared  to  the  initial  value. 
c  The  absolute  value  of  the  variation  is  given. 


TABLE  V  :  Range  of  the  <Tor2>  Dihedral  Angles  (°)  during  the  Molecular 
_ Dynamics  (MD)  Simulation _ 


Mesogen  Type 

Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

{C} 

69 

86 

36 

{C} 

62 

90 

68 

{B} 

60 

78 

117 

{B} 

66 

169 

67 

{  B } 

49 

90 

96 

The  flexibility  of  the  alkoxybenzoate  group  was  defined  as  Tori  (C-C-O-iphenyl)).  These 

dihedral  angles  (one/mesogen)  describe  the  ability  of  the  mesogens  to  rotate  with  respect  to  the 
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leader  group.  In  order  to  be  consistent  with  the  proposed  model  of  interdigitation13,16,  the 
mesogens  would  in  some  instances  need  to  be  splayed  (kinked)  with  respect  to  the  director  in  order 
to  maintain  continuous  lamellar  packing.  A  primary  location  for  this  kinking  to  take  place  may  be 
at  this  ester  bond.  Tables  VI  and  VII  indicate  that  the  disk  model  has  the  greatest  range  (=  64°, 
125°,  and  55°  for  the  cone,  disk,  and  cylinder  models,  respectively)  and  the  lowest  average 
difference  in  the  mean  dihedral  angles  (=  48°,  31°,  and  42°  for  the  cone,  disk,  and  cylinder  models, 
respectively).  The  average  range  for  the  disk  model  is  almost  twice  that  of  the  cone  and  cylinder 
models,  which  may  foretell  a  greater  tendency  towards  interdigitation  for  this  model.  It  is 
interesting  to  point  out  that  with  the  exception  of  the  disk  model  (which  has  one  gauche  bond),  all 
of  the  mean  <Tor3>  values  derived  from  the  dynamics  simulation  assume  nearly  a  transorientation 
of  approximately  180°. 

Changes  in  the  Tori,  Tor2,  and  Tor 3  torsional  angles  indicate  that  the  disk  model  exhibits 
a  larger  conformational  flexibility  than  the  cylinder  and  cone  models.  Especially  indicative  are  the 
large  range,  yet  the  relatively  small  difference  between  the  initial  and  mean  MD  dihedral  angles  for 
the  disk  model.  The  advantage  in  flexibility  for  the  disk  may  be  important  with  regard  to  the  LC 
phase  behavior  observed  for  these  materials.  In  general,  the  amount  of  interdigitation  is  primarily 

<T|Q 

governed  by  the  flexibility  of  the  leader  groups  and  of  the  attached  mesogens  .  Also  important  is 
the  persistence  of  the  interdigitation,  which  may  also  be  influenced  by  the  conformational 
flexibility.  Based  on  these  two  criteria,  the  data  may  support  the  disk  model  over  the  other  two  as 
a  more  favorable  choice  for  a  model  of  a  highly  interdigitated,  splayed  mesogenic  LC  phase. 
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TABLE  VI  :  <Tor3>a  and  <ATor3>b  Mean  Dihedral  Angles  (°)  from  the  Molecular 
Dynamics  (MD)  Simulation  in  Comparison  to  the  Initial  Values 

Mesogen 

Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

Type 

<Tor3>  <£T or3  > 

<Tor3  >  <£T or3> 

<Tor3> 

<ATor3> 

{C}  Initial 

177 

177 

179 

MD 

173  4 

184  7 

188 

9 

{ C }  Initial 

185 

170 

179 

MD 

148  37 

185  15 

177 

2 

{ B }  Initial 

119 

183 

76 

MD 

190  71 

294  111 

179 

103 

{ B }  Initial 

171 

173 

85 

MD 

202  31 

161  12 

169 

84 

{ B  }  Initial 

91 

175 

183 

MD 

186  95 

165  10 

171 

12 

a  Tor 3  =  (C-C-0-{  phenyl  1). 

b  The  mean  variation  of  the  torsion  as  compared  to  the  initial  value. 

TABLE  VII  : 

Range  of  the  <Tor3>  Dihedral  Angles  (°)  during  the  Molecular 
Dynamics  (MD)  Simulation 

Mesogen  Type 

Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

{C} 

55 

179 

48 

{C> 

103 

137 

43 

{ B} 

58 

137 

72 

{ B } 

52 

94 

41 

{B> 

52 

78 

71 

The  rotatable  bonds  between  the  two  phenyl  rings  in  the  three  {B  jmesogens  (Tori  )  as 
shown  in  Tables  VIII  and  IX,  were  also  examined.  The  results  show  little  difference  between  the 
models.  Figures  10(a)-10(c)  depict  the  variations  in  one  of  the  Tori  angles  during  the  simulation, 
showing  the  behavior  of  the  three  models  to  be  similar.  Note  that  although  the  disk,  cone,  and 
cylinder  structures  have  changed  considerably  during  the  molecular  dynamics  simulation,  they  are 
still  termed  as  such  in  order  to  be  able  to  distinguish  between  different  starting  models. 
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TABLE  VIII  :  <Tor4> *  and  <ATor4>b  Mean  Dihedral  Angles  (°)  from  the 
Molecular  Dynamics  (MD)  Simulation  in  Comparison  to  the  Initial  Values 


Mesogen 

Type 

Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

<Tor4> 

<JSFor4> 

<Tor4> 

<£Tor4> 

<Tor4> 

<ATor4> 

{B} 

Initial 

0 

0 

-2 

MD 

9 

9 

6 

6 

6 

8 

{B> 

Initial 

0 

0 

10 

MD 

-6 

6 

-12 

12 

-2 

12 

{B} 

Initial 

0 

0 

-6 

MD 

-4 

4 

-9 

9 

-12 

6 

aRotatabIe  bond  between  the  two  phenyl  rings  in  the  three  (B)mesogens. 
b  The  mean  variation  of  the  torsion  as  compared  to  the  initial  value. 


TABLE  IX  :  Range  of  the  <Tor4>  Dihedral  Angles  (°)  during  the  Molecular 
_ Dynamics  (MD)  Simulation _ 


Mesogen  Type 

Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

{B} 

46 

43 

35 

{B} 

78 

88 

72 

{B> 

87 

91 

73 

The  final  set  of  dihedral  angles  examined  were  those  within  the  siloxane  ring  (Tor5  =  ( Si- 
O-Si-O(l)  and  0-Si-0-Si(2}).  Tables  X  and  XI  indicate  the  large  amount  of  flexibility  and 
conformational  variation  present  in  the  siloxane  ring.  Once  again,  the  disk  model  exhibits  the 
greatest  flexibility,  further  supporting  the  assumption  that  the  disk  may  be  the  most  favorable 
conformation  in  a  LC  phase. 
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TABLE  X  <77>r5>a  and  <ATor5>b  Mean  Dihedral  Angles  (°)  from  the  Molecular 

Dynamics  (MD)  Simulation  in  Comparison  to  the  Initial  Values 
"Type  Isolated  Isolated  Isolated 

_ Cone  Disk  Cylinder  _ 


Value 

Eg&i-ssa 

Value 

Difference 

Value 

Difference 

I 

Initial 

60 

■w 

MD 

182 

122 

292 

2 

287 

100 

2 

Initial 

59 

76 

177 

MD 

183 

124 

54 

22 

52 

125 

1 

Initial 

299 

299 

285 

MD 

66 

127c 

64 

125° 

181 

104 

2 

Initial 

299 

59 

64 

MD 

292 

7 

297 

122c 

175 

III 

1 

Initial 

300 

293 

58 

MD 

70 

130c 

182 

111 

57 

1 

2 

Initial 

60 

66 

297 

MD 

304 

116c 

184 

118 

310 

13 

1 

Initial 

297 

297 

208 

MD 

173 

124 

315 

18 

178 

30 

2 

Initial 

64 

59 

173 

MD 

176 

112 

61 

2 

170 

3 

1 

Inital 

300 

68 

298 

MD 

64 

124c 

179 

111 

276 

22 

2 

Initial 

161 

295 

56 

MD 

302 

41 

173 

122 

143 

87 

3  Siloxane  ring  torsions  :  ( Si-0-Si-0Type  1  and  0-Si-0-SiTypc  2- 
b  The  mean  variation  of  the  torsion  as  compared  to  the  initial  value. 
c  The  absolute  value  of  the  variation  is  given. 

TABLE  XI  :  Range  of  the  <Tor5>  Dihedral  Angles  (°)  during  the  Molecul;  r 
_ Dynamics  (MD)  Simulation _ 


Isolated  Cone 

Isolated  Disk 

Isolated  Cylinder 

53 

69 

84 

56 

79 

111 

58 

68 

73 

62 

84 

64 

63 

113 

43 

64 

78 

39 

69 

60 

76 

65 

59 

72 

49 

77 

91 

47 

89 

94 

27 
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3.2.  Dynamics  on  Molecular  Pairs 

Molecular  dynamics  simulations  were  also  carried  out  on  molecular  pairs.  Intermolecular 
ordering  patterns  considered  for  the  global  molecular  topologies  were  based  on  the  “interdigitation 
of  mesogens"  structural  supposition  based  on  experimental  data13.  In  particular,  the  cylinder  and 
the  disk  models  were  examined  at  this  stage,  since  the  cone  model  was  shown  to  be  less  stable  than 
the  disk  and  cylinder. 


3.2.1  Cylinder  Model 

Previous  X-ray  diffraction  measurements  made  on  these  cyclic  siloxane-based  LC's 
indicate  a  primary  d-spacing  of  -  25  A13.  This  distance  apparently  corresponds  to  the  partial 
interdigitation  of  the  (B)  and  (C)  mesogens  of  the  proposed  cylinder  model  (cf.  Figure  11). 
These  molecules  mark  the  starting  point  of  the  molecular  dynamics  simulation.  Results  of  such  a 
simulation  may  indicate  to  what  degree  the  (B)  and  (C)  mesogens  remain  interdigitated,  if  at  all. 
Two  different  simulations  were  performed,  discussed  in  detail  below. 


3.2.1. 1  CASE  I  -  Unconstrained  Molecules 

The  starting  structure  for  this  simulation  is  the  model  shown  in  Figure  1 1.  Two  cylinder 
molecules,  each  minimized  to  remove  large  distortions  in  geometry  and  bad  steric  interactions, 
were  moved  next  to  one  another  and  their  mrermolecular  energies  were  minimized.  Less  than  five 
van  der  Waals  close  contacts  (between  2.6  A  and  3.0  A)  were  found  in  the  final  interdigitated 
structure.  A  molecular  dynamics  simulation  (details  are  presented  in  Table  lb)  was  performed, 
with  the  resulting  lowest  energy  and  statistically  averaged  structures  from  the  dynamics  run  shown 
in  Figure  12.  Because  no  constraints  were  placed  on  the  siloxane  rings  in  these  molecules,  they 
do  not  remain  parallel  to  one  another.  This  freedom  of  rotation  might  not  be  as  energetically 
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Figure  1 2:  Model  structure  of  the  unconstrained  cylinder  interdigitated  pairs  (CASE  I)  after  a  10 
pico  second  molecular  dynamics  simulation.  The  upper  structure  is  the  lowest  energy  relationship 
adopted  by  the  molecules  during  the  simulation.  The  lower  structure  is  the  statistically  averaged 
structure  of  the  molecule  pairs  during  the  simulation. 
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favorable  in  the  LC  phase  due  to  constraints  from  other  nearby  molecules  in  the  system.  This  will 
be  dealt  with  by  constraining  the  siloxane  rings  (CASE  II  of  this  section).  Nonetheless,  it  is  most 
interesting  to  note  the  large  degree  of  interdigitation  among  the  biphenyl  and  cholesteryl  mesogens. 
The  extent  of  the  interdigitation  in  this  structure  (after  a  10  pico  second  dynamics  simulation)  is 
comparable  to  that  seen  in  the  minimum  energy  starting  structure.  This  result  lends  support  to  the 
model  which  was  proposed  by  Banning  et  al.13,  suggesting  that  the  pendant  mesogens  may 
indeed  be  intetdigitated  in  the  LC  phase,  thus  accounting  for  the  primary  d-spacing  observed  with 
X-ray  diffraction  studies.  Although  a  short  simulation  time  was  used,  the  total  dynamics  run 
(including  equilibranon)  was  140  ps,  showing  no  tendency  of  the  system  to  dissociate. 

Statistical  data  for  the  interior  mesogens  torsion  angles  was  also  examined.  In  general,  it  is 
shown  that  the  torsional  angle  Tori  (Table  XII)  varies  in  a  range  of  30°  to  60°  for  the  interior 
mesogens  of  this  twin  cylinder  model.  Figures  13(a)  and  13(b)  depict  an  example  of  the  variation 
in  one  of  the  Tori  angles  during  the  simulation  (respectively  for  each  one  of  the  imerdigitated 
molecules).  Surprisingly,  there  is  not  a  remarkable  difference  in  the  conformational  flexibility  or 
stability  of  the  imerdigitated  mesogens  when  compared  to  the  results  for  an  isolated  cylinder.  In 
both  cases,  the  differences  in  the  dihedral  angles  are  small  and  the  range  similar,  except  for  one  of 
the  torsions  which  changes  from  an  initial  value  of  69'  to  a  mean  value  of  181°  during  the 
simulation.  The  average  range  for  the  imerdigitated  mesogens  is  slightly  lower  due  to  possible 
interactions  introduced  upon  packing. 

The  statistical  data  for  Tor3  shown  in  Table  XID  reveals  little  difference  for  this  dihedral 
angle  between  imerdigitated  structure  and  that  of  the  isolated  cylinder,  and  only  two  of  the  angles 
are  changing  markedly  during  the  simulation.  In  both  cases  there  is  a  tendency  of  these  angles  to 
assume  values  of  ca.  180°. 
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_  (b) 

(  )  (b)  correspond  to  the  two  mteniigitated  cylind 
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Table  XII:  Mean  Dihedral  Angles  and  Range  of  Variations  for  the  <Torl>a 
Torsion  for  the  Unconstrained  Cylinder  Model 

_ Isolated  Cylinder 

Mesogen  Initial  Mean  of  Difference  Range  Difference  Range 
Tyne  Dynamics 


(C) 

185 

177 

8 

39 

9 

41 

{C} 

300 

292 

O 

u 

36 

17 

56 

{B} 

69 

181 

112 

40 

6 

38 

{B} 

61 

55 

6 

46 

3 

84 

{B} 

60 

56 

4 

59 

2 

58 

a  See  footnote  a  of  Table  0. 


Table  XIII:  Mean  Dihedral  Angles  and  Range  for  the  <Tor3>a  Torsion  for  the 
Unconstrained  Cylinder  Model 

Isolated  Cylinder 


Mesogen 

Type 

Initial 

Mean  of 

Dynamics 

Difference 

Range 

Difference 

Range 

{C} 

179 

193 

14 

344 

116 

41 

(C> 

179 

181 

2 

43 

12 

71 

{B} 

85 

197 

112 

70 

103 

72 

f  B } 

183 

187 

4 

55 

9 

48 

{B} 

76 

258 

182 

36 

2 

43 

a  See  footnote  a  of  Table  VI. 


The  range  during  the  MD  simulation  for  Tori  and  Tort  is  shown  in  Figures  14(a)  and 
14(b)  for  the  left  molecule,  and  similarly  in  Figures  15(a)  and  15(b)  for  the  right  molecule  of  the 
interdigitated  system.  Interestingly,  the  ranges  of  change  in  Tort  are  relatively  small  in  both 
molecules  (Figures  14(b)  and  15(b)),  showing  no  large  change  in  the  amount  of  interdigitation 
These  ranges  are  again  very  similar  to  those  exhibited  by  the  isolated  cylinder  model. 

On  the  other  hand,  larger  changes  are  observed  in  Tor5  for  the  cylinder  pair  than  for  the 
isolated  cylinder  model  This  is  apparent  when  comparing  the  <ATor5>  for  example  for  the  eight 
first  angles  shown  in  Table  XIV  to  the  corresponding  values  shown  in  Table  X  and  XI,  with  the 


average  magnitudes  being  larger  by  35°  and  18°  for  the  left  and  right  siloxane  rings,  respectively. 
The  average  range  is  larger  by  27°  for  the  right  siloxane  ring  than  that  for  the  isolated  cylinder  case. 
This  indicates  that  the  siloxane  ring  undergoes  a  large  amount  of  torsional  changes  during  the 
simulation  when  the  mesogens  are  intendigitated. 

Table  XIV:  Mean  Dihedral  Angles  and  Range  of  Variations  for  the  <TorS>a 
Torsions  for  the  Unconstrained  Cylinder  Model 


Initial 

Dynamics 

Left 

Mean 

Difference 

Left 

Dynamics 

Right 

Mean 

Difference 

Right 

Range 

Left 

Range 

Right 

187 

181 

6 

176 

11 

61 

59 

177 

181 

16 

193 

16 

67 

61 

285 

335 

50 

325 

40 

92 

133 

64 

317 

253 

346 

282 

78 

127 

58 

147 

89 

147 

89 

109 

97 

297 

279 

18 

291 

6 

88 

60 

208 

281 

73 

261 

53 

52 

148 

173 

59 

114 

45 

128 

46 

130 

a  See  footnote  a  of  Table  X. 


3.2.1.2  CASE  II  -  Constrained  Molecules 

The  starting  structure  of  the  model  used  in  this  simulation  is  identical  to  the  starting 
structure  used  in  CASE  I  except  for  the  presence  of  atom  constraints.  To  eliminate  the  ability  of 
the  siloxane  rings  to  change  their  relative  orientation,  atom  constraints  were  applied  to  all  atoms 
making  up  the  siloxane  rings.  These  constraints  forced  the  selected  atoms  to  remain  fixed  at  their 
initial  minimized  {x,y,z)  coordinates.  Thus,  by  constraining  the  rings  it  was  expected  that  the 
effects  of  an  LC  phase  would  be  better  modeled.  The  details  of  the  molecular  dynamics  calculation 
are  presented  in  Table  Ic. 
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The  low  energy  and  the  statistically  averaged  structures  are  shown  in  Figure  16.  In  each 
case  the  degree  of  interdigitation  is  less  than  in  CASE  I,  as  the  pendant  mesogens,  especially  { B  J , 
move  outward  from  the  central  region  between  the  siloxane  rings.  Because  the  degree  of 
interdigitation  in  CASE  H  is  less  than  in  CASE  I,  another  dynamics  simulation  is  to  be  performed 
in  which  two  interdigitated  molecules  are  surrounded  by  additional  molecules  in  order  to  create 
boundary  conditions  for  the  inner  two  molecules  which  would  better  describe  the  LC  state. 

An  examination  of  the  changes  in  the  torsions  during  the  simulation  reveal  differences 
between  CASE  I  and  CASE  U  of  the  dynamics  simulation.  Table  XV  shows  the  statistical  data 
from  the  interior  mesogens  for  the  constrained  cylinder  model.  The  average  range  of  flexibility  for 
this  model  of  ca.  90°  is  much  larger  than  either  the  isolated  or  unconstrained  cylinder  models  This 
model  also  has  a  larger  average  difference  between  initial  and  mean  dihedral  angles  (=  61°). 
Figures  17(a)  and  17(b)  depict  the  variation  in  one  of  the  Tori  angles  during  the  simulation  (for 
each  one  of  the  interdigitated  molecules)  in  which  a  large  change  was  observed,  clearly  seen  also  in 
Figure  16.  By  constraining  the  siloxane  ring,  larger  torsional  movements  of  the  leader  group  near 
the  siloxane  ring  (Tori)  were  observed,  while  the  range  and  average  difference  for  the  Tor3 
torsions  were  very  similar  to  those  of  the  unconstrained  case  because  it  was  several  atoms  away. 
Contrary  to  this  behavior,  the  Tor2  values  did  not  change  as  much  as  in  the  unconstrained  case 
because  of  the  hinderance  to  ring  movement.  The  Tor4  and  TorS  values  are  very  similar  to  the 
unconstrained  and  isolated  cylinder  model  values  due  to  the  siloxane  rings  being  hindered. 
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figure  17:  Statistical  data  of  the  10  pico  second  molecular  dynamics  simulation :  trace  of  the 
variation  in  a  dihedral  angle  Tori :  (a)  and  (b)  correspond  to  the  two  i.itetdigitated  cylinders. 
These  angles  correspond  to  those  shown  in  Figure  13. 
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Table  XV:  Mean  Dihedral  Angles  and  Range  of  Variations  for  the  <Torl>a 
Torsion  for  the  Constrained  Cylinder  Model 
_ Isolated  Cylinder 


Mesogen 

Type 

Initial 

Mean 

Dynamics 

Difference 

Range 

Difference 

Range 

(Cl 

185 

82 

103 

129 

9 

41 

{C} 

300 

300 

0 

39 

17 

56 

IB) 

69 

79 

10 

38 

6 

38 

<B} 

61 

115 

54 

187 

3 

84 

{B} 

60 

200 

140 

44 

2 

58 

a  See  footnote  a  of  Table  II. 


3.2.2  Disk  Model 

3.2.2.1  CASE  III  -  Constrained  Molecules 

As  in  the  cylinder  model,  in  which  partial  interdigitation  of  the  mesogens  corresponds  to 
the  primary  d-spacing  observed  in  X-ray  diffraction  studies,  the  partial  interdigitation  of  these 
mesogens  in  the  disk  model  may  also  correspond  to  the  observed  X-ray  data.  Molecular  dynamics 
calculations  were  performed  on  two  interdigitated  disk  molecules,  each  minimized  to  its  minimum 
energy  conformation  as  seen  in  Figure  1 8.  As  was  done  before,  the  inrermolecular  energy  between 
the  two  molecules  was  minimized  and  less  than  five  van  der  Waals  close  contacts  were  found  in  the 
final  structure.  Atom  constraints  were  also  imposed  m  this  simulation  in  order  to  maintain  the 
initial  relationship  between  the  two  siloxane  rings.  All  atoms  in  the  siloxane  rings  were 
constrained  to  their  initial  minimized  (x,y,z)  coordinates.  The  details  of  this  simulation  are  given 
in  Table  Id.  The  low  energy  and  the  statistically  averaged  structure  of  this  system  after  a  10  pico 
second  simulation  are  shown  in  Figure  19. 

The  range  during  the  molecular  dynamics  simulation  for  Tor2  and  Tor5  is  very  small  since 
the  rings  are  constrained.  Statistical  data  for  the  Tori  torsions  show  relatively  large  variations 
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Figure  19:  Model  structure  of  the  constrained  disk  intendigitated  pain;  (CASE  DI)  after  a  10  pico 
second  molecular  dynamics  simulation.  The  upper  structure  is  the  lowest  energy  relationship 
adopted  by  the  molecules  during  the  simulation.  The  lower  structure  is  the  statistically  averaged 
structure  of  the  molecule  pairs  during  the  simulation. 
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analogous  to  the  constrained  cylinder  model.  Indeed,  it  is  interesting  to  note  that  the  so  called  disk 
model  transforms  into  a  conformation  more  resembling  the  cone  model  of  Figure  2.  A  more 
extensive  dynamics  simulation,  similar  to  the  one  mentioned  in  the  previous  section  should  be 
conducted  here  in  order  to  provide  additional  boundary  conditions  to  the  molecules.  Also,  a  pair  of 
unconstrained  disk  models  should  be  studied  in  order  to  examine  if  the  siloxane  rings  absorb 
energy  in  a  manner  analogous  to  the  unconstrained  cylinder  model.  Finally,  it  should  be  noted  that 
the  disk  model  proposed16  is  somewhat  different  than  the  disk  model  examined  here.  The 
proposed  model  has  a  unit  director  associated  with  it  as  the  mesogens  have  been  splayed  so  that 
they  point  in  one  direction.  The  disk  system  examined  here  is  more  analogous  to  a  radial  system 
with  the  mesogens  radiating  in  all  directions 
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3.3.  X-ray  Diffraction  Simulations 

Simulated  X-ray  diffraction  patterns  of  the  two  cylinder  pair  systems  (CASES  I  and  II)  and 
of  the  disk  pair  system  (CASE  III)  were  examined.  The  simulated  meridional  scattering  patterns 
for  the  unconstrained  cylinder  model  (CASE  I)  are  shown  in  Figures  20(a)-20(d).  Additional 
reflections  appear  for  the  final  structure  as  compared  to  the  initial  structure.  The  origin  of  the  large 
peak  in  these  scans  at  28=1.405°,  d=62.8A  may  be  due  to  the  relatively  small  numbers  of  units 
interdigitated.  It  has  been  shown  that  longer  units  have  an  effect  on  the  resulting  X-ray  pattern16 
due  to  effects  of  the  noninterdigitated  end  mesogens.  Scattering  patterns  for  the  constrained 
cylinder  model  (CASE  II)  are  shown  in  Figures  21(a)-21(d).  Once  again,  there  is  considerable 
change  in  the  scattering  pattern  between  the  first  and  last  structure  of  the  dynamics  simulation. 
Table  XVI  summarizes  the  d-spacings  of  these  structures  for  CASES  I  and  II.  Clearly,  an  increase 
in  the  number  of  peaks  in  the  meridional  scans  indicates  that  CASE  II  is  a  more  ordered  structure. 
The  equatorial  scattering  patterns  for  the  disk  pair  model  (CASE  III)  are  shown  in  Figure  22. 

Table  XVI:  X-ray  Diffraction  Maxima  (A)  from  Calculated  Meridional  Scattering 
Patterns  of  the  First  and  Last  Structures  (CASES  I  &  II) 

CASE  I  CASE  II  Experimental13 


First 

Last 

First 

Last 

62.8 

62.8 

62.8 

62.8 

15.7 

20.9 

20.9 

20.9 

24.0 

10.4 

12.6 

12.6 

11.9 

7.0 

7.9 

7.9 

8.0 

5.7 

6.3 

6.0 

4.8 
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Figure  20:  CASE  I:  Meridional  scattering  section  for  (a)  the  starting  structure  of  the  10  pico 
second  MD  simulation,  (b)  the  final  structure  of  the  10  pico  second  simulation,  (c)  the  lowest 
energy  structure  of  the  10  pico  second  simulation,  (d)  the  statistically  averaged  structure  of  the  10 
pico  second  simulation.  Incident  X-ray  beam  is  perpendicular  to  the  siloxane  rings. 
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(a) 


(c) 


Figure  21:  £&aEJL  Meridional  scattering  section  for  (a)  the  starring  structure  of  the  10  pico 
second  MD  simulation,  (b)  the  final  structure  of  the  10  pico  second  simulation,  (c)  the  lowest 
energy  structure  of  the  ten  ps  simulation,  (d)  the  statistically  averaged  structure  of  the  10  pico 
second  simulation.  Incident  X-ray  beam  is  perpendicular  to  the  siloxane  rings. 
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Figure  22:  CASE  III  :  Equatorial  scattering  section  for  (a)  the  starting  structure  of  the  10  pico 
second  MD  simulation,  (b)  the  final  structure  of  the  10  pico  second  simulation  (c)  the  lowest 
energy  structure  of  the  tlO  pico  second  simulation,  (d)  the  statistically  averaged  structure  of  the  10 
pico  second  simulation.  Incident  X-ray  beam  is  perpend  to  the  siloxane  rings. 
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Section  4 

CONCLUSIONS 


Molecular  dynamics  simulations  offer  new  insights  into  the  conformational  flexibility  of 
these  cyclic  siloxane-based  LCs.  Interdigitarion  between  the  cholesteryl  4'  -  aiiyloxybenzoate  and 
biphenyl  4'  -  aiiyloxybenzoate  mesogens  pendant  on  the  siloxane  ring  is  certainly  present  in  the 
simulated  structure,  but  a  quantitative  measure  of  the  interdigitarion  is  still  to  be  calculated, 
possibly  by  the  application  of  boundary  conditions  to  the  system.  The  isolated  disk  model 
exhibited  the  most  flexibility  and  greatest  stability  as  indicated  by  mean  dihedral  angles  and  range 
for  certain  key  torsions.  In  general,  all  three  isolated  models  displayed  large  conformational 
flexibility,  which  may  be  important  with  regard  to  the  LC  phase  behavior.  Results  from  the 
dynamics  simulation  of  the  cylinder  pairs  indicate  the  large  conformational  flexibility  of  the 
siloxane  rings.  Movements  of  the  interdigitated  mesogens  were  much  higher  for  a  fixed  rings 
system  than  for  the  case  where  the  ring  dihedral  angles  were  free  to  move.  X-ray  scattering  pattern 
calculations  for  the  structures  generated  during  the  dynamics  run  demonstrate  higher  order  than  for 
the  initial  models,  as  shown  in  Table  XVI. 
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